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ABSTRACT 

We designed a simulation of an intracellular 
differential equation model of the dynamics of 
malaria with immune control and treatment which 
considered malaria parasites in the liver and blood. 
We considered transmission dynamics of malaria and 
the interaction between the infection in the liver and 
blood. The disease free equilibrium of our model was 
asymptotically stable when the basic reproduction 
number is less than one and unstable when it is greater 
than one. Numerical simulations show that if the 
immune response is strong with effective treatment, 
malaria infection will be cleared from an infectious 
human host. A treatment strategy using highly 
effective drugs against malaria parasites with strong 
immune response can reduce malaria progression and 
control the disease. 

Keyword: Mathematical model, Malaria parasite, 
Hepatocyte, Erythrocyte, Meroziote, Sporozoite, 
Immune response, Treatment. 

I. INTRODUCTION 

Malaria is a life threatening mosquito borne blood 
disease caused by a plasmodium parasite and children 
are particularly susceptible to the disease. In 2015, an 
estimated 306,000 children under 5 years of age were 
kills mostly in the African region (WHO World 
Malaria Report, 2015). Once transmitted to the human 
by a blood feeding Anopheles mosquito, the parasites 
initially multiply in the human liver, before they 
progress to the pathologic blood stage. Immediately 
the parasite (sporozoites) first enters the human host, 
there is a pre- erythrocytic development. After 
inoculation into a human by female anopheles 
mosquito, sporozoites invade hepatocytes in the host 
liver and multiply there for 5-12 days, forming 
hepatic schizonts. These then burst, liberating 
merozoites into the bloodstream where they 


subsequently invade red blood cells. These blood 
infections can last for months, and only once sexual 
precursor cells, the gametocytes have matured, the 
malaria parasite are able to leave the human host and 
to continue the life cycle in the insect vector. In the 
mosquito midgut, the parasite are able to differentiate 
into their sexual forms, the female macrogametes and 
male microgametes, and to then undergo sexual 
reproduction in order to newly combine their 
chromosomal sets. The midgut phase lasts for 
approximately 20 hours and includes two phases of 
stage conversion, the rapid conversion gametocyte 
into fertile gametes upon activation and the 
conversion of zygotes into motile and invasive 
ookinates that once formed, immediately exit the gut 
lumen by traversing the midgut epithelial cell layer. 
Subsequently, the ookinates settle down at the basal 
site of the midgut epithelium and convert to sessil 
oocysts in which sporogonic replication takes place. 
This replication phase requires roughly 2 weeks and 
results in the formation of infective sporozoites that 
migrate to the salivary glands to be released into the 
human dermis with the next bit of the mosquito 
wherewith the life cycle of plasmodium is completed 
(Aly et al, 2009; Ghosh and Jacobs-Lorena, 2009; 
Kuehn and Pradel, 2010; Menard et al, 2013; Bennink 
et al, 2016). 

Sexual precursor cells the intraerythrocytic 
gametocytes develop in the human blood in response 
to the stress factor (Pradel, 2007;Kuehn and Pradel, 
2010). A time period of about 10 days is required for 
gametocyte development in P. falciparum, during 
which they pass five morphological stages. Once the 
gametocytes mature and is ingested with the blood 
meal of an Anepheles mosquito, they are activated in 
the mosquito midgut by environmental stimuli, and 
gametogenesis is initiated. Signals inducing gamete 
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formation include a drop of temperature by 
approximately 5°C which is mandatory for 
gametocyte activation and the presence of the 
mosquito derived molecule Xanthurenic acid (XA), a 
metabolic intermediate of the tryptophan catabolism. 
An additional trigger of game to genesis is the 
increase of extracellular pH from 7.2 to about 8 
(Kawamoto et al, 1991;Billker et al, 1997; Garcia et 
al, 1998; Sologub et al, 2011). 

The periodic bouts of fever that occur in the malaria 
coincide with the synchronized rupture of 
plasmodium-infected red blood cells. This causes the 
release of parasites en masses into the blood stream, 
along with pigments and toxins that have accumulated 
inside the red blood cells as a result of the parasites 
metabolic activities. The presence of large quantities 
of parasite material in the blood triggers a dramatic 
immune response, mediated by the secretion of 
cytokine modules by the cells of the immune system 
(Hommel and Gilles, 1998). Some cytokines such as 
tumor necrosis factor (TNF), interferon gamma, 
interleukin 12 and interleukin 18 enhances the immune 
response, stimulating macrophages and other immune 
cells to destroy parasites by phagoytosis and by the 
production of toxins. Other cytokines include 
interleukin 4, interleukin 10 and TGF-beta help to 
regulate the immune response by dampening these 
effects (Malaguarnera and Musumeci, 2002). 

II. Related Literatures 

Chi-Johnston (2012) develops and analyze a 
comprehensive simulation model of P. falciparum 
within-host malaria infection and transmission in 
immunologically-naive humans. There model 
incorporates the entire lifecycle of P. falciparum 
starting with the asexual blood stage forms 
responsible for disease, the onset of symptoms, the 
development and maturation of sexual stage forms 
(gametocytes) that are transmissible to Anopheles 
mosquitoes, and human to mosquito infectivity. The 
model components were parameterized from malaria 
therapy data and from arrange of other studies to 
simulate individual infections such that the ensemble 
is statistically consistent with the full range of patient 
responses to infection. Human infectivity was 
modeled over the course of untreated infections and 
the effects were examined in relation to transmission 
intensity expressed in terms of the basic reproduction 
number. Adamu (2014) developed a mathematical 
model to study the dynamics of malaria disease in a 
population and consideration were given to the 


interaction between the parasites and the host (human 
beings), such that the susceptible and the infected 
classes were allowed to interact freely without 
quarantining any of the either classes. In their model, 
first order equation that describes the dynamics of the 
susceptible class and the infected class under the 
influence of the parasite was used. The result of the 
qualitative and stability analysis showed that if 
preventive measure is not put in place, the susceptible 
and infected classes will reach a stable equilibrium 
point which can be disastrous to the population and 
they recommended specific measures of controlling 
the disease. 

Johansson and Leander (2010) used three 
compartment of susceptible, infectious and recovered 
in their work and they showed that the recovered are 
neither quarantined nor removed from the entire 
population rather they enter the susceptible class 
again. Tabo et al (2017) developed a mathematical 
model which considers the dynamics of P. falciparum 
malaria from the liver to the blood in the human host 
and then to the mosquito. There results indicated that 
the infection rate of merozoites, the rate of sexual 
reproduction in gametocytes, burst size of both 
hepatocytes and erythrocytes are more sensitive 
parameters for the onset of the disease. They 
suggested that a treatment strategy using highly 
effective drugs against such parameters can reduce on 
malaria progression and control the disease. There 
numerical simulation shows that drugs with efficacy 
above 90% boost healthy cells and clear parasites in 
human host. However, all these models are limited to 
treatment, non considered treatment and immune 
response. Here, we formulated a more detailed model 
to study the intracellular dynamics of malaria with 
immune control and treatment using mathematical 
model. Our aim is to study the interaction between 
malaria and immune response with treatment measure 
through a mathematical model and carry out a 
sensitivity analysis to determine the parameters that 
controls the disease. 

m. MODEL FORMULATION 

3.1. THE ASSUMPTIONS OF THE MODEL 

The disease is spread by transmission through 
mosquito to human interaction; 

By immunology memory, the immunity of infected 
/infectious individuals might be rapidly restored when 
they are re exposed to the infection. 
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Individual can loss their immunity when they are not susceptible. 

continuously exposed to the parasite and go back to Treatment can either be successful or fail. 



3.2. Variables (Compartments) 

The model is made up of ten (10) compartments 
which comprises of (x), Uninfected hepatocytes (liver 
cells), (p),Free sporozoites (malaria parasites in the 
liver), (y), Infected hepatocytes, (Ty), Treated 
infected hepatocytes, (/?), Recovered hepatocytes, 
(Pi), Free merozites (malaria parasite in the blood), 
(5), Uninfected erythrocytes (red blood cell), (/), 
Infecfed erythrocytes, (T,), Treated infected 
erythrocytes and (R^, Recovered erythrocytes. 

Parameters 

xp recruitment level of uninfectec hepatocytes 
a-L natural death rate of both uninfected, infected, 
and recovered hepatocytes 
/? rate at which hepatocytes are being infected 
p death rate of malaria parasites (sporozoites) 


a 3 

a 7 


a 


a A 


a 5 


a. 


0 


rate at which free sporozoite is inoculated into 
the hepatocyte by mosquitoes 
treatment rate of infected hepatocytes 
movement rate of treated hepatocytes to 
recovered class 

natural death rate of erythrocytes (red blood 
cells) 

recruitment level of erythrocytes from bone 
marrow 

rate at which the uninfected erythrocytes are 
being infected 

rate at which the infected erythrocytes produce 

free parasites (merozoite) 

disease induced death rate of infected 

erythrocytes 

disease induced death rate of infected 

hepatocytes 

the rate at which infected hepatoctesproliferate 
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9 the rate at which infected erythrocyte 

proliferate 

d 2 rate treatment of the infected erythrocytes 

d 3 movement level of infected erythrocytes to the 
recovered class 

ky rate at which infected hepatocyte produces 
meroziotes 

b 3 recovered red blood cells due to immune 

response 

b A recovered liver cells due to immune response 

b movement rate of the recovered hepatocytes to 
susceptible class. 

b x movement rate of the recovered red blood 
cells to susceptible class 

3.3. The model equation 

dx 

— ip — aqx — /?xp + bR 


dt 


dy 

dt 


dp 

dt 


= a 3 y - pp 


= fixp + 4>y - $y - a i y ~ °-ey ~ kypi ~ b 4 y 


dT„ 


dB 

dt 


O 7 tii Ty diTy 

+ b 4 y (3.1) 


dt ~ - 1 y 1 y 

dR 

—— — diT v — ciiR — bR 
dt 1 y 1 


— Y — aBp! — a 2 B + 


dPi 

dt 


- a 4 I — pp! 


dl 

— = aBp! + 91 — d 2 I — a 2 I — a 5 I — b 3 I 
dT, 

~dt ~ d 2 I ~ dsTl ~ azTl 


dR t 

dt 


— d 3 T[ cl 2 Ri biRi + b 3 I 


Let the initial conditions be 

x(0) = x o ,y(0) = y 0 ,T y ( 0) = T y0 ,R( 0) = R 0 ,B( 0) 


negatively much faster than that of the cell 
concentration. 

Notice from equation (3.1) that the production rate of 
the parasite (p), from the livercells is proportional to 
the rate at which they are removed and are at 
equilibrium, i.e., a 3 y — pp — 0. So we let 

a 3 y 

P =- 

p 

Also from equation (3.1), we observe that the rate of 
production of the parasite(p 1 ), from the red blood 
cells is proportional to the rate at which the are 
removed and are at equilibrium, i.e., a 4 / — pp± — 0. 
So we let 

CL^I 

a Vl=T a, 

Substituting p = ^-and Pi = ~ into equation (3.1) 

reduces the model to eight non linear ordinary 

differential equations and this will make the 

quantitative analysis much easier. Now we rewrite the 

equations as: 

dx a 3 y 

— — ip — a ± x — Bx -1- bR 

dt P 

dy a 3 y a 4 I 

-^ = P x ~ + (f>y - Z y - a i y - a-sy - ky — 

- b 4 y 


dT v 


(y tl \ Ty d ] Ty 


dt ^ ~ 1 ^ 

dR 

—— — diT v — a t R — bR 
dt 1 y 1 

+ byy (3.3) 

dB claI 

—— — y — aB - a 2 B + biRi 

dt P 

dl aj 

— — aB -h 91 — d 2 I — a 2 I — a s I — b 3 I 

dt p 

dT, 

— d 2 1 — d 3 T, — a 2 Ti 

dRi 

~ d 3 T[ — ti 2 Ri — b 1 R 1 + b 3 I 


= B o ,I(0) = I o ,T I (0) = T IO ,R 1 (0) 

— Rio (3-2) 

3.4. Equilibrium state analysis 

The equilibrium state is the uninfected state and for 
malaria infection to manifest, the individual must be 
bitten by an infected mosquito. Also, the rate of 
change in sporozoites and merozoites concentration 
will be positively much faster than that of the cell 
concentration and for it to clear, the rate of change in 
sporozoites and merozoites concentration will be 


Because the model s are items of populations and in 
two interacting cell population, that is, the liver cells 
which produces sporozoites and the red blood cells 
which produces merozoites. The liver cell and the red 
blood cell population size at time t are respectively 
represented as 

x(t) + y(t) + T y (t ) + R(t) 

— N{t)andB(t ) + /(t) + T;(t) 

+ R 1 {t)=N 1 (t) 
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3.5. Existence and Positivity of solutions 


At this point we let the time derivative of 


Having that all the parameters in equation (3.1) are A L (t)and A r (t) along solutions of system (3.2) for 


non negative, we assume a stable population with per 
capita recruitment of susceptible liver cells, 
susceptible red blood cells, death rate of liver cells 
both natural and disease induced, death rate of red 
blood cells both natural and disease induced. At this 
point we normalize the population size of both the 
liver cells and red blood cells to one (1) each and 
show that the system is epidemiologically and 
mathematically well-posed in the feasible region T 
given by 


where 

Ai 


T = A, xA,c 


x 


and A r 


= {(*,y,T y )6R3:N<^} 

= f(B, [,T,j e R+: N, < j-J 


Theorem 1: There exists a domain T in which the 
solution set [x,y, T y ,B, I, T,] is contained and 
bounded. 

Proof: Given the solution set [x, y, T y ,B, I, Tj] with 
positive initial conditions (3.2), we let the liver 
population be represented as 

x + y + Ty + R — 1 (3.4) 

=> R = 1 — x — y — T y 

while the red blood cell population is represented as 
B + / + T[ + R 4 

= 1 (3.5) 

=>R 1 = 1-B-I-T I 

Omitting the equation for RandR 1 in our analysis 
gives equation (3) as 

dx _ a 3 y 


— = xb — a.iX — Bx- 
dt r 1 B 


dy a 3 y 

—— — fix -+ <py ~ ty ~ a 4 y 

dt b 


+ b(l -B -I -7)) 

<* 6 y 


a 4 I 
ky — 
B 


dT v 


-b 4 y 


dt 

= <y- a x T y 

~ djy 

dB claI 

— — y — aB - a 2 B + b 4 { 1 — B — I — 7)) 

dt b 

dl cla.1 

— — aB - V 61 — d 2 I — a 2 I — a s I — b 3 I 

dt b 

dT, 

~dt ~ d 2 I ~ dsTl ~ a{Yl 


liver cells and red blood cells respectively be 
calculated thus, 

A L (t) = x(t) + y(t) 

+ T y (t ) (3.7) 

a 3 y 

A L (t) = ip — a x x - fix -1- b( 1 — B — I — 7)) 


+ x 

+0y - <y - %y - 


%y 

b 

a 6 y 


a a I 


where 


ky- + (y- ai T y 
d i T y - b 4 y 


A l = x + y + T y 


Remember that in the absence of the 
diseas ed-^Ty^y , <py, b 4 yanda 6 y will be equal to 
zero. Then we obtain 

A L (t) = xp — a 4 x — a 4 y — a 4 T y + b{ 1 — A L ) 
A L (t) = xp ~ %(x + y + T y ) + b(l - A L ) 

A L (t) = xp- a 1 A L + b - bA L 
AlW + (a i + b) A L 

+ b (3.8) 

We shall integrate both sides of equation (3.8) using 
integrating factor method according to (Kar and Jana, 
2013; Birkhoff and Roffa, 1989) to obtain: 

A l + P{t)dt — F(t) 

A l < e ~f p ^ dt (| e^ p(t)dt F(t)dt + c) 

where P(t ) — a 4 + b and F(t) — xp + b. Let the 
integrating factor be 

r(t) = e J p ( f )dt — g/Oi +b)dt _ g (a 1 +fc)t 


e (a 1 +fc)t 


Then integrating equation (3.8) by inputting r(t) = 
gives 

/lift) <A-(J r(t).F(t)dt + C ) 

(3.6) =, At ( t) < (| » + b)dt + c) 

1 


1 


A L (t) < r M e(ai+b)t + C 

LK J gCai+Ot \ ( a i + b) 
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OP + b) 


< 


(% + b ) 

+ Ce~^ +h)t 


(3.9) 


Here, C is the constant of integration and if we let t 
oowe have that 

(xp + b ) 

AM = (^Trr x+y+T ? 

But 

x < — 

Also, 

A r (t) 

= B(t) + I(t) + Tj(t) 

A r (t) = y- aB — - a 2 B + b 4 ( 1 -B-I-Tj ) 

H' 

a A 
+ aB — 

B 

+91 — d 2 I — a 2 I — a 5 I — b 3 I + d 2 I — d 3 Tj — a 2 Tj 


(3.10) 


(3.11) 


where 


A r — B + I + Ti 


Also, in the absence of the disease, 
91, a 5 I, b 3 Iandd 3 Tj will be zero. Then we have 
A r (t) = y~cl 2 B - a 2 I - a 2 T, + b t ( 1 - A r ) 
A r (t) —y — a 2 {B + / + 7)) + b 4 (l — A r ) 

A r (t) = y - a 2 A l + b 4 - b 4 A L 
A r (t) + (a 2 + h 1 )A r 

<Y + b 1 (3.12) 

Using integrating factor method on equation (3.12), 
we have 
A r (t) 

(y + b i) 


< 


(a 2 + h) 

+ C 1 e“ (a2+i,l)t 


(3.13) 


Here, C 

is the constant of integration and if we let t -> cowe 
have that 


A r (t) = 


(y + h i) 


— B + I + T, 


is positively invariant and the systems for the liver 
cells and red blood cells are respectively well-posed 
epidemically and mathematically. Then for the initial 
starting point A L G R+ andd r G R+ the trajectory lies 
on T. Thus, we focus our attention only on the region 
T. 

3.6. Disease Free Equilibrium point 

To study the equilibrium state and analyze the 
stability of the system, we set the right side of 

equation (3.3) to zero. Thus, we have 

a 3 v 

xp — a ± x — f3x -1- bR — 0 

ct 3 y CL4I 

fix -1 - <py — ~ a i y ~ a 6y ~ ky - b 4 y 

B 6 

= 0 

Cy - a i T y - d iTy = 0 
d-J’y — a x R — bR + b 4 y 
= 0 (3.15) 


CL a I 


aB —— a 2 B + b^Ri = 0 
B 


aB + 91 — d 2 l — a 2 l — a s I — b 3 I — 0 
B 

d 2 I — d 3 T, — a 2 T\ — 0 
d 3 T 1 — a 2 R 1 — b 1 R 1 + b 3 I — 0 

If we label equation (3.15) as (3.15i) to (3.15viii), 
then (3.15ii) gives 


/?x 


a 3 y 


a a. I 


(% + b 4 ) 

But 

Y 

B < — 
a 2 

Observe from the dynamics describe by the systems From (3.15i) we get 
(3.2), (3.10) and (3.14) that the region 


+ 0y - <y - a i y - a +y - ky - b 4 y = o 

B B 

( a 3 a 4 I \ 

f/?x-1- (p — ( — a 4 — a 6 — k - b 4 )y — 0 

=> y — 0 

From (3.15iii) we have 

Cy “ a i T y ~ d i T y = 0 
But y — 0, then we have 

(oq + dP)Ty — 0 Ty — 0 

From (3.15iv) we obtain 

d x T y — a t R — bR + b 4 R — 0 
Since T y — y — 0, we have 

(% + b)R = 0 => R — 0 
From (3.15i) we have 

a 3 y 

xp — a 4 x — /3x -1- bR — 0 


(3.14) 


But y and R are all equal to zero, then we have 

xp 

xp — ayx = 0 => x = — 


CL -1 


a a. I 


T = f(x,y,T y , B, I, T 

,) G E+: N < — 

: N! < —| 


a i 

a 2 ) 


aB -1 - 91 — d 2 I — a 2 I — a 5 I — b 3 I — 0 

B 
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&4 


>) /= 


aB -1 - 9 — d 2 — a 2 — a 5 — b 3 

[A 

=> / = o 

From (3.15vii) we have 

d 2 l — d 3 T l — a. 2 T\ — 0 

But / = 0, then 

(d 3 + a 2 )Ti = 0 => T { — 0 
Substituting 7) = / = 0 in(3.15viii) we obtain 
(a 2 + b 1 )R 1 = 0 => — 0 


0 


Also, if we substitute I — R ± — 0 into (3.15v) we get 

y — a 2 B — 0 => B — — 

cl 2 

There, the disease free equilibrium point of the model 
is given as 

O = (x, y, T y , R, B, I, T,, R^) 

= (—,0,0,0, —,0,0,o) (3.16) 

3% a 2 / 


3.7. Existence and stability analysis of disease free equilibrium 

To find the Jacobian matrix of the model system, we differentiate equation (3.3) with respect to 
x, y, T y , R, B, I, Tj, R ± respectively to obtain. 

dx* 


dt 


dy* 

dt 


a l P 

x* + 

a 3 

—/?x —— 

L p J 


L ii\ 


y* + [b]R* 


\o a *y] 


p — 
L P J 

X* + 


dB* 

dt 


dR* 

dt 


do 


dT* 

dt 


I 

B 




r a 4^i 

b 4 

y* + 

-k — 


1 y 


- [fly* + [~( a i + d)]Ty 


[b 4 ]y* + [d ± ]T; + [—(a, + b)]R* 


dl* 

dt 


CL^I 

-a -a 2 

B 


B* + 


a A B 


-a- 


I* + [bi]Rl 


aj 


\ a 4 B 

a — 

B* + 

a -h 6 — d 2 — a 2 — a 3 — b 3 

L p J 


L 11 


dR{ 

dt 


dr; 

dt 


= [d 2 ]l* + [-(d 3 + a 2 )]r; 


= [d 3 ]T,* + [b 3 ]I* + [-(a 2 + MR* 


We examine the stability of the disease free equilibrium using equation (3.16) 

Pv a , 


J(Q) := 


t-ia. 


0 




|ra 


+ 4 -C-a 1 -a g -b 4 


l l + d l) 


0 


0 0 


0 0 


“( a l + b ) 


0 -a 2 

0 0 

0 0 

0 0 


aya 


|.ia„ 


0 

0 

0 

0 

-aya 4 
^a 


+ 0 - d 2 - a 2 - a 5 - b 3 


"( a 2 + d s ) 

d. 


-( a 2 + b l). 
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J(Q)-M 


-X -a ( 

0 (|) - A. - C+ ■ 

0 
0 


Pf a 

|ia 

Pfa 


3 


3 


Ha. 


■a. - a, - b . 
1 6 4 


-X - a, - d, 


-X - b - a, 0 


-X - a„ 


aya , 


0 0-X + ' 

0 
0 


Ha- 


0 

0 

0 

Ha, 


■ a 2- a 5- b 3- d 2 


4- a 2-d 3 


-A, - a -b 

2 ly 


We need to show that all the eigen values of the matrix J(Q) have negative real part. Observe that the first and 
fifth columns contain only the diagonal terms and this forms the two negative eigen values X x — — a 4 and X 2 — 
—a 2 , the other six eiginevalues can be obtained from the sub-matrix, ] 2 {Q ), formed by excluding the first and 
fifth rows and columns of J(Q). thus, we have 

f P^a 3 

<p — A. — +-- - - a & - 0 0 0 0 0 


JjlQ) := 


|ra 




b - aj - X 


cxya 


0 0 -X + 

0 

0 


l-ia. 


a 2- a 5 “ b 3- d 2 


-*-a 2 -d ; 


-X - a - b 
2 ly 


In the same way, the third and sixth column of J\(Q) contains only the diagonal term 
which forms negative eigenvalues^ = — (a 1 + b — h 4 ) and X 4 — —(a 2 + b ± ). The 
remaining four eigenvalues are obtained from the sub - matrix 


J 2 (Q) := 


-x-i ; + 


Pv a 3 

iraj 

? 

0 

0 


“ a i - a 6“ b z 


-X - - d j 


ocya / 


0 -X + 


Ha„ 


a 2"a 5 -b 3 -d 2 


-X - a - d 
2 3 J 
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The eigenvalues of the matrix J 2 (Q ) are the roots of the characteristic equation 


/ pxpa 3 \ ( aya 4 \ 

( —A -I- h (p — ( — a-, — a 6 — b 4 ) (—A — a-, — d 4 ) ( —A H-1 - 9 — d 2 — a 2 — a s — b 3 ) (—A — 

V [IQ) \ (IQ, 2 / 

- a 2 ) = 0 


which translates to 


As — — d-3 cl 2 ,A 7 — 


(hpa^ 
pa 4 


+ cp — ( — a 1 


CLc 


b 4 , 


aya 4 

and Ao — -h 6 — d 2 — a 2 — a 5 — b 3 

pa 2 


This implies that the eigenvalues A 12 j6 are both less that zero i.e., A x < 0 ,A 2 < 0, A 6 < O.If + 0 < 


^ + a 4 -|- cig and 


aya 4 

t*a 2 


MO i 


+ 6 < d 2 + a 2 + a 5 + b 3 , clearly, A 7 and A 8 will respectively be less than zero 


(A 7 < 0 and A 8 < 0) and that means that the steady state is asymptotically stable. But if 


a A and 


aya 4 
Ita 2 




+ 0 > + a 4 + 


+ 9 > d 2 + a 2 + a 5 + b 3 , A 7 and A 8 will respectively be greater than zero (A 7 > 0 and A 8 > 0), 


we conclude that the steady state is unstable. 

3.8. Basic Reproduction Number R 0 

The basic reproduction number R 0 is the average number of secondary infectious infected by an infective 
individuals during the whole cause of disease in the case that all members of the population are susceptible 
(Zhien et al, 2009; Olaniyi and Obabiyi, 2013). 


To obtain R 0 for model equation (3) we use the next generation technique (Van den Driessche and Watmough, 
2002; Diekmann et al, 1990). We shall start with those equations of the model that describes the production of 
new infections and change in state among infected liver cells and red blood cells. 


Let H = [ x, y, T y , B, I, T,\ where T denotes transpose. 


dH 

■ 1J7 = F(H) - V(H) 

\Ma 3 y 

+ 0y 


F(H) = 


\ia 4 

0 

aya 4 I 


pa 2 


+ 91 


0 


,V(H) = 


(3.17) 

kya 4 l 

0y + +-+ a 6 y + b 4 y 

F 

_ <y + ajy + djy 

d 2 I + a 2 I + a 3 I + b 3 I 
—d 2 I -I- d 3 Tj + a 2 T, 


Finding the derivatives of F and V at the disease free equilibrium point O gives F and V respectively as 


F := 


Pw, 


ira 


+ 00 


a Y a 




0 

0 

+ 00 

0 J 


V := 


(^ + a i + a 6 + b z 

< 

0 

0 


a l + d l 
0 

0 


o o ) 

o o 


a 2 + a 5 + b 3 + d 2 0 


2 3 J 
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v -> 


C+ a i + a 6 + b 4 


C-a, + C-d +a-a. + a-b. + a-d, + a -d. + b-d. + (a,l^ a l + d 1 
^ 1 -1 16 14 11 61 4 1 \ 1/ 

0 0 


a 2 +a 5 + b 3 + d 2 


( a 2 + d 3 )-( a 2 + a 5 + b 3 + d 2) a 2 + d 3 


F-V 


Pv|/a q 


4> + 


|oa 


^ + a i + a 6 + b 4 


aya 


0 + 


Ha„ 


a 2 + a 5 + b 3 + d 2 


0 J 


| F.V~ 1 \ = 0 


Here, we can obtain the basic reproduction number R 0 from the tree and determinant of the matrix FV 1 — G. 
R 0 — W (G) - - trace(G) + ^trace(G ) 2 — 4 det(G) (3.18) 


Observe that det(G) = 0, so we have 

Pxpci 3 0 


Rn = 


+ 


+ 


yaa 4 


+ 


9 


a 1( u(<( + a ± + a 6 + b 4 ) ( + a x + a 6 + b 4 a 2j u(a 2 + a 5 + b 3 + d 2 ) a 2 + a s + b 3 + d 2 


(3.19) 


From equation (3.19). 


Wa-3 
a if* 


+ <p is the multiplication ability of the disease in the liver cells and the probability 


that an individual will move to the second level which is the infection of the red blood cells; 


is the 


<^ + CL'i + CZ5 + £>4 

average duration of infectious period of the liver cells before the release of the parasite to invade red blood 


cells; 


yaa 4 

a 2 H 


+ 9 is the multiplication ability of the disease in the red blood cells and the probability that the 


individual will be infectious; 


a 2+ a S+t ) 3+d2 

Let the basic reproduction number R 0 be written as 


is the average duration of the infectious period of the red blood cells. 


where 


Rq — Rl + Rr 


R, = 


(3.20) 


Pxpa 3 


+ 


0 


Q-i /r(£ + a 4 + ci 6 + h 4 ) £ + a 4 + n 6 + b 4 

vaa 4 9 

+ 


and R r — -— , , . , , 

a 2/ u(a 2 + a 3 b 3 + d 2 ) a 2 + a 3 + b 3 + d 2 
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We have R L describing the number of liver cells that one infectious liver cell infects over its expected infectious 
period in a completely susceptible liver cell population. While R r describes the number of red blood cells 
infected by one infectious red blood cell during the period of infectiousness in a completely susceptible red 
blood cell population. 

3.9. Existence of Endemic Equilibrium point 

Endemic equilibrium point describes the point at which the disease cannot totally be eradicated from the 
population. We shall show that the formulated model system (3.3) has an endemic point and we let O** be the 
endemic equilibrium point. 


Theorem 2: the intracellular malaria model system (3.3) has no endemic equilibrium when R 0 < 1 but has a 
unique endemic equilibrium when R 0 > 1. 

Proof: Let = (x**,y**,Ty*,R**,B**,I**,Tjl f *,R{*') be a nontrivial equilibrium of the model system (3.3); 
i.e., all components of ®** are positive. If we solve equation (3.3) simultaneously having in mind that ®** =£ 0 
we have that 

d- a x R** - bR** + b 4 y** = 0 
d x T;* + b 4 y** 


_ 


** _ 

ly ~ 




a ± + b 

(y** - aj** - dj** = 0 


cl i + d 4 

Therefore (3.21) can be rewritten as 


1 (d^y** \ 

R** — -- , +b 4 y**) 

CLi + b Vaq + d 4 J 


d?I** - a 2 T',** - d?T,** = 0 


_ 

I — 


d 2 r 

cl 2 T d-i 


3 1 1 


(3.21) 

(3.22) 

(3.23) 

(3.24) 


d 2 d 3 I** 


T~) ** _ 

Hi ~ 


ccB* 


a 2 3" d 3 

\d 2 d 3 + b 3 (a 2 + ^ 3 )]/ 

(a 2 + ^i)( a 2 3 " d 3 ) 
aj** 


+ b 3 l** = (a 2 + bi)R{* 


(3.25) 


B 


+ or* - d 2 r - a 2 r - a 5 r - b 3 r = o 


CLa 


a B** -h 6 — d 2 — a 2 — a 5 — b 3 — 0 

(d ? + a ? + a q + b 3 — 9)u 

B** = — -------— (3.26) 


CC CLa 


y-a\ 


\d 2 + a 2 + a 5 + b 3 — 0 )/A a 4 I* 


CC CLa 


cl 2 


\d 2 + a 2 + a 5 + b 3 - 9 V 


CC CLa 


[d 2 d 3 + b 3 (a 2 + d 3 )]I _ 

1 (a 2 + b t )(a 2 + d 3 ) 


b![d 2 d 3 + b 3 (a 2 + d 3 )]/** - (d 2 + a 2 + a 5 + b 3 - 0)(a 2 + b^a -2 + d 3 )I* 
f (d 2 + a 2 + a 5 + b 3 - 9)^i\ 
a 2 -—- - Y 


CCCLa 


(a 2 + b 1 )(a 2 + d 3 ) 
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j** _ 


(d 2 + a 2 + a 5 + b 3 - 9)fi \ 

2 l rv n . I • 


CC CL a 


(a 2 + h 1 )(a 2 + d 3 ) 


b 1 [d 2 d 3 + b 3 (a 2 + d 3 )] - (d 2 + a 2 + a 5 + b 3 - 0)(a 2 + &i)(a 2 + d 3 ) 


(3.27) 


a 3 y 

/?x + (py** - (y** - a 4 y** - a 6 y* 


a 4 l** 


b 4 y** = 0 


a 3 a 4 /** 

/3x** - \- cp — ( — a 4 — a 6 — k v -b 4 = 0 

q k 


/ a 4 r* \ n 

X** — I <( + a 4 -I- ci 6 + ky -h b 4 — cp j —— 

V U / p CL 


(3.28) 


xp — a 4 x** — f3x* 


a 3 y 


+ bR** =0 


/ a 4 I** \ n / a 4 r* \ 

xp — a 4 ^ + a 4 + cl 6 + ky -h b 4 — cp J — - (3 ^ + a 4 + ci 6 + ky -1 - b 4 — cp J 


k 

b / ddy** \ 

H- — (-——h b 4 ) = 0 

a 4 + b \a 1 + d t J 

( a 4 I** \ 

(3 + a 4 + a 6 + ky -h b 4 — <pj 


k a 3 y 

' (3 a 3 q 


dd 


k a 3 y 


/ 3a 3 ii a 4 + b Va 4 + d 


-+i, 4 ) 


y 


— \xp — a 4 K + a 4 + a 6 + k y 


ad** 


+ h/i 


*)& 


y = 


xp - a 1 (? > + a 4 + a 6 + k y —\- b 4 — xp'j 


pa 3 


+ a 4 + cig + k 


1X4,1 

y~ 


) b 1 

f d 

) a 4 + b 1 

\a 4 + d t 


(3.29) 


3.10. Sensitivity Analysis of the Basic Reproduction Number R 0 

Observe that the basic reproduction number R 0 is in the form R 0 — R L + R r , where R L and R r are functions of 
nine parameters respectively. But R 0 is a function of sixteen parameters which comprises of the basic 
reproduction number at the liver site and the basic reproduction number at the blood site. To control the 
disease, these parameter values must control R 0 , such that its value will be less than one (R 0 < 1). Therefore 
change in the parameter values, results in change in R 0 and if we let 

<3l ~ (fit a l> a 3> a 6’ k> 0) an d (Jr ~ (y> a > a 2 ' a 4> a 5< k> b 3> d 2 , 0) 


then the rate of change of R 0 for a change in the value of parameter q can be estimated from a normalized 
sensitivity index 


7^0 _ 

L q ~ 


dR L Rl 
dq L ' (Jl 


dR r R r 
dq r q r 


(3.30) 



dR L Rl r 
-.— and Z„ r 

Hr 


dq L kh 


dR r R r 

dq r q r 
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Table 1 Parameter values for calculating R 0 and Numerical Simulation of the Model 


Parameters 

Description 

Value 

Range 

Reference 

0 

Recruitment levelof uninfected hepatocytes 

3 x 10 8 

Mota et al, 2001 

a 4 

Natural death rate of both uninfected, infected, 
treated and recovered hepatocytes 

0.002- 

0.0067 

Mota et al, 2001 

P 

Rate at which hepatocytes are infected 

4 x 10~ 9 

Tabo et al, 2017 

R 

Natural death rate of malaria parasite 

48 

Tabo et al, 2017 

a 3 

Rate at which infected hepatocytes produce free 
sporozoites 

0.181 

Esteva et al, 2009 

b 

Rate at which recovered hepatocytes move to 
susceptible class 

1.3 

x 10~ 4 

Chitnis, 2008; Mohammed and 
Orukpe, 2014 

< 

Treatment rate of infectious hepatocytes 

0.95 

Mohammed and Orukpe, 2014; 
Castillo-Riquelme et al, 2008. 


Movement rate of treated infectious hepatocytes 
to recovered class 

0.1 

Ducrot et al, 2008 

a 2 

Natural death rate of erythrocytes (Red blood 
cells) 

0.0083 

Anderson et al, 1989 

y 

Recruitment level of erythrocytes 

2.5 x 10 8 

Austin et al, 1998 

a 

Rate at which uninfected erythrocytes are being 
infected 

2 x 10 10 

Dondorp et al, 2000 

a 4 

Rate at which infected erythrocytes produce free 
merozoites 

16 

Chiyaka et al, 2010 

a 5 

Disease induced death rate of infected 
erythrocytes 

0.24 

Chiyaka et al, 2010 

a 6 

Disease induced death rate of infected 
hepatocytes 

2.0 

Tabo et al, 2017 

0 

The rate at which infected hepatocytes 
proliferate 

3 x 10“ 5 

Estimated 

9 

The rate at which infected erythrocytes 
proliferate 

2.5 

x 10~ 5 

Estimated 

d 2 

Rate at which infectederythrocytesare being 
treated 

0.95 

Mohammed and Orukpe, 2014 

d 3 

Movement rate of treated infected erythrocytes 
to recovered class 

0.01 

Ducrot et al, 2008 

ky 

Rate at which infected hepatocytes produce 
merozoites (malaria parasite) 

16 

Tabo et al, 2017; Chiyaka et al, 
2010 

^3 

Recovered erythrocytes due to immune response 

4.56 

Estimated 

b 4 

Recovered hepatocytes due to immune response 

0.0035 

Shah and Gupta, 2013 

b 1 

Movement rate of recovered erythrocytes to 
susceptible class class 

1.37 
x 10~ 4 

Chitnis, 2008; Molineaux and 
Gramiccia, 1980 


To calculate the value of R 0 , we use the parameters as stated in table 1. 

— Rl + R r 

D _ /?0% 0 

k l - 


+ ■ 


Rr - 


a-ifitf + cti T cig + b 4 ) £ + cl 4 + n 6 + b 4 

yaa 4 9 


+ 


a 2 /u(a 2 +a 5 + b 3 + d 2 ) a 2 + a 5 + b 3 + d 2 
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4 x 10~ 9 x 3 x 10 8 x 0.181 3 x 10~ 5 


R, = 


+ 


0.004 x 48(0.95 + 0.004 + 2 + 0.0035) 0.95 + 0.004 + 2 + 0.0035 

= 0.382512257 


2.5 x 10 8 x 2 x 10“ 10 x 16 


Rr = 


10 


+ 


2.5 x 10" 5 


r 0.0083 x 48(0.0083 + 0.24 + 4.56 + 0.95) ' 0.0083 + 0.24 + 4.56 + 0.95 

= 0.348723952 

R 0 = 0.382512257 + 0.348723952 = 0.731236209 

The normalized sensitivity index of the basic reproduction number with respect to /?, xp, a 1( a 3 , 
a 6 , q, 0 is given by 

Rl Ml Rl 

qi dq L ' q L 


= ^L^L = ( _ ^P a 3 

13 dp ' p + a 4 + a 6 + b 4 )J V P 


)(!)-■ 


0110964 x 10 15 


jRl _ Rl^ _ ( _ P a 3 _\ / R_L 

^ dxp ' xp \a 4 q(( + a 4 + a 6 + b 4 )J \xp 


■K 


69 x 10~ 18 


z: l = 


dR, R 


“ 3 da , 


■-=(- 

\CL 


pxp 


Rl _ dR L R l _ 
------ 


q(Z + a 4 + a 6 + b 4 )J \a 3 
1 \ /R 




4018101352 


d<p (p \( + a 4 + a 6 -I- b 4 / V cp 


m)- 


4311.2116885 

0 

Jfll 3a 4 ' a 4 V( a iM^ + a iR + + a-^qb^ 2 ' (<" + a x + a 6 + h 4 ) 2 y 

= -9156.8606658 


_ 


Pxpa 3 tf + 2a 4 + a 6 + b 4 ) 


+ 


l) 




Rl _ 3R L _ 
------ 


dR L 

Rl _ ( 

( pxpa 3 \ 

dq 

■ R ~ 1 

Ka-^q 2 ^ + a 4 + a 6 + b 4 )J 


d( ( 


pxpa 3 


■ + ■ 


A^\ 

— = - 0 . 

R / 

, )(,) 


0488762883 
R, 


r l = djh Rl 
1,4 db 4 ' b 4 


Z « L = 9 Rl Rl 
° 6 3a 6 ' a 6 


a iM(( + a 4 + a 6 + h 4 ) 2 (<( + a 4 + a 6 + b 4 ) 2 J V £ 

= -0.0521115791 

pxpa 3 (p \(R l 


+ 




%/!«■ + a 4 + a 6 + h 4 ) 2 (( + a 4 + a 6 + b 4 ) 2 J \b 4 

= -14.144571463 

pxpa 3 (p \(R l 


+ 


OS 


a 4 q(y + a 4 + a 6 + b 4 ) 2 (<( + a 4 + a 6 + b 4 ) 2 J \a 6 

= -0.0247530001 
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The normalized sensitivity index of the basic reproduction number with respect to y, a, a 2 ,a 4 ,a 5 ,p, b 3 ,d : 
given by 


_ dRr Rr 
Ql dq r ' q r 


^R r _ dR r R r _/ aa 4 

Y dy ' y \a 2 n(a 2 +a 5 +b 3 + d 2 )J \y 




96 x 10- 18 


7 R r _ dRr Rr _ 


ya 4 


) = 3.0401721 x 10 18 

a4 da 4 'a 4 Va 2 p(a 2 + a 5 + b 3 + d 2 )) (a 4 ) 


da ' a \a 2 ii(a 2 + a 5 + b 3 + d 2 ) 

^ Rr dR r R r ( ya \ /R 


.0004750269 


7 R r _ dR r Rr _ 

6 ~ dd ' 6 


)©- 


jR l _ dRr Rr _ 

“ 2 da 2 ' a 2 


+ 


2422.4090582 

9 


jR r _ dR r R r _ 


(a 2 + a 5 + b 3 + d 2 ) 

yaa 4 (2a 2 + a 5 + b 3 + d 2 ) 

(ajfi + a 2 a s n + a 2 b 3 y. + a 2 d 2J u) 2 ' (a 2 + a 5 + b 3 + d 2 ) 2 / 

= -1761.6460452 

yaa 4 


Rr\ 

a 7 ) 


dfi n 


jR r _ dR r Rr _ 
as da 5 ' a 5 


a 2 ju 2 (a 2 + a s + b 3 + d 2 ), 
yaa 4 


)© — 


0000527808 


9 


+ ■ 


m 


jR r _ dRr Rr 
1,3 db 7 ' b 7 


-( 


a 2 p(a 2 +a 5 + b 3 + d 2 ) 2 (a 2 + a 5 + b 3 + d 2 ) 2 / Va 5 

= -0.0879950061 
yaa 4 


9 


+ ■ 


OS 


_ 


dR r R r 


d dn d,'} 


(a 


a 2 ^(a 2 + a 5 + b 3 + d 2 ) 2 (a 2 + a 5 + b 3 + d 2 ) 2 / Vb 3 

= -0.0046313161 

yaa 4 9 \ /R 


2 d ( a 2 + a 5 + b 3 + d 2 ) 2 + (a 2 + a 5 + b 3 + d 2 ) 2 7 Vd 2 
= -0.00222303174 


Table 2: The effect of the parameters on R L . 


Parameters 

Value Range 

Effect on R 01 

P 

4 x 10~ 9 

9.0110964 x 10 15 

0 

3 x 10 8 

1.69 x 10“ 18 

a 3 

0.181 

4.4018101352 

0 

3 x 10- 5 

4311.2116885 

CL 

0.004 

-9156.8606658 

R 

48 

-0.0488762883 

( 

0.95 

-0.0521115791 

b 4 

0.0035 

-14.144571463 

a 6 

2.0 

-0.0247530001 


, 9 is 
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Table 3: The effect of the parameters on R r . 


Parameters 

Value Range 

Effect on R 01 

y 

2.5 x 10 8 

1.69 x 10“ 18 

a 

2 x 10“ 10 

3.0401721 x 10 18 

a 4 

16 

0.0004750269 

e 

2.5 x 10" 5 

2422.4090582 

a 2 

0.0083 

-1761.6460452 

\x 

48 

-0.0000527808 

a 3 

0.24 

-0.0879950061 

b?, 

4.56 

-0.0046313161 

d 2 

0.95 

-0.0222303174 


The sensitivity index Z((3), Z(ip), Z(a 3 ) and Z(<p) are all positive and this shows that the value of R L increases 
as the value of p,xp,a 3 and <p increases. The remaining indices Z(a 1 ),Z(h 4 ),Z(a 6 ),Z( / u) and Z(() are 
negative, indicating that the value R L decreases as a 1; b 4 , a 6 , /i and ( increases. Actually, the effectiveness of 
control may be measured by its effect on R L . if the reduction in R L < 1 can be maintained by the parameters 
a lt b 4 , a 6 , fx and then it will reduce the endemicity of the disease. This implies that these parameters can help 
in reducing the rate of malaria infection over time in the liver and if it is maintained, the transmission of the 
disease may decrease, causing the cases in the liver population to go below an endemicity threshold. 

Similarly, the sensitivity index Z(y),Z(a),Z(a 4 ) and Z(0) are all positive indicating that the value of R r 
increases as the value of y, a, a 4 and 6 increases. The indices of remaining parameters 
Z(a 2 ),Z(b 3 ),Z(a 5 ),Z(ju) and Z(d 2 ) are negative, and this shows that the value of R r decreases as 
a 2 , b 3 , a 5 , ye and d 2 increases. Since the effectiveness of control may be measured by its effect on R r and if the 
reduction in R r < 1 can be maintained by the parameters a 4 , b 4 , a 6 , jx and then endemicity of the disease in 
the erythrocyte will be reduced. Therefore, the parameters a lt b 4 , a 6 , jx and ( can help in reducing the rate of 
malaria infection over time in the erythrocyte and if maintained, the transmission of the disease may decrease, 
causing the cases in the erythrocyte population to drop beyond the endemicity threshold. 

IV. Numerical Analysis and Results 

The numerical behavior of system (3.3) were studied using the parameter values given in table 1 and by 
considering initial conditions, (p = {x(0),y(0),T y (0),R(0),S(0),/(0),'r / (0),R 1 (0)}. The multiplication 

ability of meroziote in the hepatocyte is _j- 0 — 1.13128, while the probability that the red blood cell will 
be infected by sporozoites is --= 0.3381234151. Also, The multiplication ability of sporoziote in the 

erythrocyte is _|_ g — 2.0080571285, while the probability that the human host will be infectious is 

a 2 n 

---= 0.1736623656. 

a 2 +a s + b3 + d 2 

The numerical simulation are conducted using Matlab software and the results are given in figure 2-4 where 
figures 2i - 2iii illustrate the behavior of the reproductive number R L for different values of the model parameter 
b 4 and b 3 respectively. Figures 3i - 3iii also show the behavior of the reproductive number R r for different 
values of the model parameter ( and d 2 respectively where ( is represented with g 4 . Lastly, figures 4i - 4viii 
and 5i - 5viii show the varying effects of the immune system and treatment controls. 

The basic reproduction number of the system is given by 

R 0 = R l + R r = 0.731236209 < 1 
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indicating that the basic reproduction number is less than one. Therefore, the disease free equilibrium is stable 
showing that malaria infection can be controlled in the population using adequate treatment method. However, 
it also confirms the result of the sensitivity analysis of R L and R r in tables 2 and 3 respectively. We then state 
that with effective treatment of infectious human, the future number of malaria infection cases will reduce in 
the population. 





b4=strong immunity in Liver cell 

b3=strong immunity in Red Blood cell 

b4,b3 =lmmune Responses 

fig 21 

fig 2ii 

fig2iii 


Figures (2i - 2iii): Numerical Simulation of the Basic Reproduction Number R_L and R_r using different rate 

of b_4 and b_3, (Immune Control) 



g4=Treatment in Liver cell 
fig 3i 



d2=Treatment in Blood cell 
fig Ii 



g4,d2 ^Treatment 
fig3iii 


Figures (3i - 3iii): Numerical Simulation of the Basic Reproduction Number R_L and r using different rate of 

g_4 and d_2 (Treatment Control) 


@ IJTSRD I AvailableOnline@www.ijtsrd.coml Volume-2 I Issue - 6 ISep-Oct2018 


Page: 377 





















International Journal of Trend in Scientific Research and Development (IJTSRD) ISSN: 2456-6470 



t=time fig4vii t=time fig4viii 


Figures (4i - 4viii): Numerical Simulation of model system (3.3), when there are immune and treatment control 

from 0 -8days. 
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t=time fig5vii t=time fig5viii 


Figures (5i - 5viii): Numerical Simulation of model system (3.3), when there are immune and treatment control 

from 0 -30days. 


The numerical simulation of immune response, 
treatment and disease free equilibrium point were 
performed to establish long term effects. Parameter 
values used in the simulations are given in table 1. The 
simulation of the basic reproduction number, 
R_0=R_L+R_r as in figures 2i - 2iii shows that the 
immune response is effective in reducing the density 
of the parasites both in the liver cells and red blood 
cells. It indicates that the infection rate of the 
hepatocytes and erythrocytes are respectively reduced 
as the merozoites are suppressed and the sporozoites 
being cleared. Also, from figures 3i - 3iii, we observe 
that there is a perfect treatment since the reproduction 
numbers R_L and R_r under treatment are all less 


than one and R_0=R_L+R_r is less than one. This 
implies that there exists the clearance of malaria 
parasites in both the liver and blood. Therefore with 
this reduction in infectious reservoirs, malaria can be 
greatly reduced in the population. If efficacy was 
equal to zero, that is R_0=(R_L+R_r )>1, immune 
response and treatment would have been useless. 

Figures 4i - 4viii and 5i - 5viii show the disease free 
dynamics of malaria infection at hepatocytes of the 
liver and erythrocytes of the blood. The result shows 
that in the absence of malaria, the susceptible 
hepatocytes and erythrocytes respectively increases. 
Also, there is a sharp fall in the density of the 
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infectious hepatocytes and erythrocytes. This 

indicates that the population of merozoites and 
sporozoites in the hepatocytes and erythrocytes will 
respectively decrease. Observe that new infectious 
mosquitoes repeatedly bite an individual to continue 
the life cycle to naive individuals, activating the 
immune response against the infection. With an 
increase in treatment effectiveness, the density of the 
uninfected and recovered hepatocytes and 

erythrocytes increases, while the population of 

infected hepatocytes and erythrocytes decreases to 
lower value because the efficacy of the treatment used 
is high. 

V. Discussion and Conclusion 

The proposed study of the simulation of an 
intracellular differential equation model of the 

dynamics of malaria with immune control and 
treatment was designed and analyzed using ten 
compartments which were later simplified to eight 
compartments. The model studied malaria infection 
both in liver and blood. It also incorporated the effect 
of immune response and treatment of the infection 
respectively in the liver and blood stages. The 
analysis of the model as was presented by the 
positivity and existence of the systems solution shows 
that solutions exist. The results in this model indicate 
that the disease free equilibrium is asymptotically 
stable when R_0=(R_L+R_r )<1 and unstable when 
R_0=(R_L+R_r )>1. in this study, the parameters, 
Cb_3,b_4,and d_2 were significant in the successful 
clearance of malaria parasites. The sensitivity indices 
of these parameters were negative which indicates 
that increase in them results to reduction in malaria. 
The simulation result shows that with effective 
treatment, the density of uninfected hepatocytes and 
erythrocytes, treated hepatocytes and erythrocytes and 
recovered hepatocytes and erythrocytes increases. 
This simply means that the number of merozoites in 
the liver and sporozoites in the blood will be reduced 
and this implies clearance of malaria. 
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